Equation of state and collective frequencies of a trapped Fermi gas along the 

BEC-unitarity crossover 
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We show that the study of the collective oscillations in a harmonic trap provides a very sensitive 
test of the equation of state of a Fermi gas near a Feshbach resonance. Using a scaling approach, 
whose high accuracy is proven by comparison with exact hydrodynamic solutions, the frequencies 
of the lowest compressional modes are calculated at T — in terms of a dimensionless parameter 
characterizing the equation of state. The predictions for the collective frequencies, obtained from the 
equations of state of mean field BCS theory and of recent Monte-Carlo calculations, are discussed 
in detail. 
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The availability of Feshbach resonances and the con- 
sequent possibility of tuning the interatomic potential 
through the application of a magnetic field have recently 
stimulated much experimental and theoretical work on 
the BEC-BCS crossover in ultracold trapped Fermi gases 
0. Several relevant and already measured quantities, 
like the release energy, the density profiles and the fre- 
quency of the collective oscillations, are related in a sen- 
sitive way to the equation of state of the homogeneous 
system which varies significantly through the crossover. 
For small and negative values of the scattering length a 
the equation of state approaches the limit of a non in- 
teracting Fermi gas while, for small and positive values, 
bound molecules can be formed and, at zero temperature, 
the system behaves like a dilute molecular Bose-Einstein 
condensed gas. In the intermediate regime important 
theoretical issues remain. For example the intermolecu- 
lar scattering length on the BEC side exhibits a non triv- 
ial dependence on the free atom scattering length 0, [j| . 
Furthermore such interactions should give rise to beyond 
mean field corrections in the equation of state 0, not 
accounted for by the Bogoliubov approximation. Similar 
corrections could also result from the composite nature of 
the molecules. Finally, near resonance, where the scat- 
tering length becomes larger than the average distance 
between particles, no simple many-body approach is ap- 
plicable to this strongly correlated system. Numerical 
calculations of the equation of state of a uniform inter- 
acting Fermi gas at T=0 along the crossover have been 
carried out in the past along the mean field BCS (MF- 
BCS) approach and, more recently, through ab initio 
simulations based on Monte Carlo (MC) algorithms 0,0- 

Since the frequencies of the collective oscillations can 
be measured with high precision, it is of major interest 
to investigate their dependence on the equation of state 
along the crossover. From a careful and systematic anal- 
ysis of the frequencies one should be able to infer the 



actual equation of state. It has already been pointed out 
[8( that the collective frequencies of a T=0 superfluid 
Fermi gas, trapped in a harmonic potential, approach 
well defined values in the important BEC and unitarity 
limit regimes, where the density dependence of the chem- 
ical potential can be inferred from general arguments. In 
the intermediate region various investigations, based on 
the use of the hydrodynamic theory of superfluids and 
suitable parametrizations of the equation of state, have 
appeared recently [rnil2lll^ll^ [la[^ll7l|. In the mean 
time first experimental results |3, on the frequencies 
of the lowest axial and radial compression modes on ultra 
cold gases of 6 Li across the Feshbach resonance have also 
become available. 

The purpose of this paper is to establish the accuracy 
needed in the measured frequencies in order to determine 
experimentally the equation of state. With this goal in 
mind we will focus on two specific theories for the equa- 
tion of state which exhibit a rather different density de- 
pendence of the chemical potential through the crossover. 
These are the MC calculations and the MF-BCS approx- 
imation discussed above. These equations of state are 
then employed within a hydrodynamic approach in order 
to calculate the relevant collective frequencies. Hydrody- 
namic theory is justified at T=0 by the superfluid nature 
of the system and has been already successfully employed 
in trapped Bose-Einstein condensed gases. We will re- 
strict here the discussion to the case of positive scattering 
lengths, where the experimental conditions for achieving 
the superfluid hydrodynamic regime are less severe, in 
contrast with the BCS regime where the smallness of the 
gap gives rise to stringent constraints 16]. The determi- 
nation of the collective frequencies, for a given equation 
of state, is by no means a trivial task because, in gen- 
eral, the solutions must be found numerically in the non 
symmetric configurations relevant for experiments. Fur- 
thermore since the frequencies vary only moderately with 
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the scattering length, the analysis should be carried out 
with high precision. In the following we will discuss only 
the lowest relevant modes, namely the monopole mode 
for spherical geometry and the lowest axial compression 
mode for cigar configurations. 

Let us first start discussing the two equations of state 
and, in particular, the dependence on density n of the 
chemical potential (J,(n). Instead of the density, we use 
the Fermi wavevector kp defined by kp = 3ir 2 n. We 
restrict ourselves to the case of a wide Feshbach reso- 
nance where the scattering length a is the only relevant 
length associated with the two-body interaction, as it is 
indeed the case for the dominant resonance in 6 Li 
and also for 40 K. This leads to introduce the dimension- 
less coupling parameter a, which varies from zero to 
infinity when one goes from unitarity to the bosonic limit 
of molecules. Similarly we express the chemical poten- 
tial in units of the Fermi energy Ep = h 2 k F /2m. The 
equations of state are displayed in Fig. 1. 
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FIG. 1: Chemical potential fi(n) (shifted by half the molecular 
binding energy h 2 /2ma 2 ) in units of the Fermi energy Ef as a 
function of the coupling parameter l/fc^a. The Monte-Carlo 
data are extracted from Ref. [7]. 

At unitarity the result n/Ep — 0.44 extracted from 
MC calculations 0, is in fair agreement with experi- 
ments 0, |2(J, while the MF-BCS theory gives a larger 
result n/Ep = 0.59. The slopes near unitarity are very 
similar, but the approach to the BEC limit is quite dif- 
ferent reflecting the fact that MC calculations are con- 
sistent with the recent prediction om = 0.6 a for the 
molecule- molecule scattering length 3|, while the MF- 
BCS 5] yields the significantly larger value a« = 2a. 

In order to evaluate the frequencies of the collective 
modes, different routes have been followed in the litera- 
ture. One class of methods assumes that the equation of 
state of the gas can be closely approximated by a poly- 
tropic expression /i(n) ~ n 7 and makes use of the analyt- 
ical results known [21j for this case. Note that the poly- 
tropic approximation is exact in the BEC regime where 
fi = gn, as well as at unitarity where dimensionality con- 
siderations yield /i cx n 2 / 3 . Other approaches are based 
on exact numerical solutions of the hydrodynamic equa- 
tions |23l ■ Here we will consider another method based 



on the use of scaling transformations. Scaling has been 
already successfully employed to describe the dynamical 
behaviour of Bose- Einstein condensates [U |3| ■ We will 
see that the results of the scaling approach agrees remark- 
ably well with exact solutions of hydrodynamic equations 
in all cases of practical interest we have considered. The 
advantage of the proposed method is its full applicability 
to anisotropic traps where exact numerical solutions are 
not always available. On the other hand it can be also 
applied to regimes where the equation of state is not of 
the polytropic form. 

In our scaling approach we use the following Ansatz 
for the space dependence of the density distribution 

n(xi,x 2 ,X3,t) = 7(t)no(a 1 (t)x 1 ,a2(t)x2,a 3 (t)x 3 ) (1) 

where no(r) is the equilibrium density and 7(f) = 
Yli a i(t)- The factor j(t) ensures that the total parti- 
cle number Jdrn(r,t) is conserved. Correspondingly, 
we take u s ,i(r) = — (di/a^ij for the superfluid veloc- 
ity field which ensures that the equation of continuity 
n(r) + V [n(r) v s (r)] = is automatically satisfied. We 
will then insert the above Ansatz for the density and for 
the velocity field into the hydrodynamic Lagrangian 
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dr 



mn(r) v 2 (r) — e(n(r)) — n(r)V/ l (r) 



(2) 



of our zero temperature superfluid. Here m is the atomic 
mass and the chemical potential /x(n) is linked to the 
energy per unit volume e(n) by n(n) = de(n)/dn. The 
density dependent energy term e(n) in the Lagrangian 
has been here obtained in the Local Density Approxima- 
tion, consistently with the hydrodynamic description. 

In the following we will specialize to the case of small 
oscillations, suited to the collective modes in the linear 
regime. Writing <Zj = 1 + where £j -C 1, in the equa- 
tions of motion obtained from the Lagrangian @, we 
have to carry out the calculation up to first order in e^. 
After some algebra we find the following simple eigen- 
value equations for the mode frequencies 10: 



where we have naturally used e'i 



LO 2 £j 



(3) 



and we have 



defined the relevant dimensionless parameter 



r = 



dni 



2 (Vho) 
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(4) 



where V/j (r) = J2 i (l/2)moj 2 x 2 is the trapping harmonic 
potential with x\ — x, X2 = y and x% — z. The brack- 
ets are for the equilibrium average over the gas cloud, 
for example (Vho) — J dr no(r) T4 (r) and the density 
profile is evaluated in the Thomas-Fermi approximation 
n{n) + V/i (r) = (j,q, consistently with the hydrodynamic 
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description. Notice that the parameter T does not de- 
pend on the anisotropy of the trap. Using the virial the- 
orem it is easy to prove that, for a polytropic equation 
of state /u(n) oc n 7 , the parameter T coincides with the 
exponent 7 so that at unitarity one has T — 2/3, while in 
the BEC regime r = 1. In general the coefficient T is a 
function of the dimensionless combination kp (0)a where 
fcp(O) is fixed by the central value of the Thomas- Fermi 
density distribution according to k F (0) — 37r 2 n(0). The 
actual dependence of T on kp (0)a is determined only by 
the equation of state n(n). Finally the above derivation 
can be also reformulated as starting from an exact vari- 
ational principle for the hydrodynamic frequencies. 

In the isotropic case u>i = il, in addition to the 
quadrupole mode ui 2 = 2Q 2 , Eq.© gives the result ui 2 = 
(3r + 2) fl 2 for the monopole mode and, in particular, the 
well-known result tu — \fb il in the BEC limit r — 7 — 1 . 
In the axisymmetric geometry lj\ = L02 — u>± = u) z , 
we find, in addition to the quadrupole mode u> 2 — 2lu 2 _, 
the two m = compressional solutions: 

^ = \ [2(r + i)^i + (r + 2)^ 

±y/[2(T + 1) u\ (r + 2) ujlf + 8T 2 Lo\ Z 2 } (5) 

For the relevant case of cigar geometry lo z <C uj±_ one finds 
oj 2 = uj 2 (3 — l/(r + 1)) for the lowest axial compression 
mode, giving w/w 2 = \/5/2 in the BEC limit. For the 
radial compressional mode one finds instead to 2 = 2(T + 

We will now show that the scaling result is extremely 
accurate by comparing its predictions with essentially 
exact solutions of the hydrodynamic equations. Our 
method [13] makes use of known exact analytical results 
for a large class of model equation of states fi a nai(n). 
Then the actual u(n) is very closely approximated by 
one member of this class. Moreover the effect of the 
small difference between u(n) and fianaiin) is taken into 
account using a perturbative approach (this is called 
the "corrected" model). The absolute precision of this 
method has been checked j2^| to be of order 1CP 3 at 
least. The relevant ingredients needed to calculate the 
mode frequencies are the reduced gas density in the trap 
n(r) = n(r)/no and the corresponding chemical potential 
p,(r) = u(r)/uo, both normalized to their central values 
no and A class of soluble models is described by the 
law n = fi p exp Kfe_o Pkfi k ], where p and are param- 
eters. Here we consider only the cases K = (which is 
the polytropic model h = pP) and K = 1, because they 
give already by far a good enough precision. The model 
K = 1 is called " quasi-polynomial" . In this last case the 
solution is not fully analytical, but it involves only find- 
ing the proper root of a low order polynomial, which is a 
quite trivial numerical task [22j |. 

We have carried out the calculations of the collective 
frequencies in two different physical situations. The first 



one is the isotropic geometry which is particularly con- 
venient for deriving exact solutions of the hydrodynamic 
equations. The second one is the elongated cigar geom- 
etry employed in many experimental set-up. Note that 
the results v 2 — 4 for the unitary case and v 2 = 5 for the 
BEC limit are much more separated in the spherical case 
than in the cigar geometry where one finds instead, for 
the lowest axial mode, v 2 =2.4 and v 2 = 2.5 respectively. 




atan (1/k p (0) a) 



FIG. 2: Reduced square v 2 — uj 2 /ft 2 of the lowest monopole 
frequency in isotropic trap for the Monte-Carlo (MC) and 
the mean-field BCS (MF-BCS) equation of states as a func- 
tion of l/&F(0)a, where fcf-(O) is the Fermi wavevector at the 
trap center. Full lines corresponds to the results of all our 
methods, which can not be distinguished at this scale. In- 
sert : (uncorrected) quasi-polynomial (dotted-dashed line), 
corrected polytropic (full line) and scaling (dashed), all com- 
pared to the corrected quasi-polynomial model. 

We display our results in Fig. 2 for the spherical ge- 
ometry. We have considered the two equations of state 
provided by MC calculations and by the MF-BCS ap- 
proximation. The accuracy of the scaling approach is in 
both cases remarkable since at the scale of this figure 
it coincides with the exact result. The excellent quality 
of the scaling approach can be understood because it is 
exact not only for the polytropic case, but also if one 
includes first order corrections to the polytropic equa- 
tion of state [27| , as can be proven from the variational 
formulation mentioned above. In order to display the 
precision of our results, we have plotted the results in 
the insert of Fig. 2 on a much larger scale, taking the 
MC equation of state as an example. More precisely 
we have taken as a reference the frequency v 2 pcorr ob- 
tained by the corrected quasi-polynomial method, which 
we believe is the most accurate. Indeed we see first that 
the corrected polytropic model and the corrected quasi- 
polynomial are in remarkable agreement, since the differ- 
ence between them is at most 8.10 -4 . Hence the conver- 
gence of our method toward the exact result is almost al- 
ready achieved at the level of the simpler corrected poly- 
tropic model. This clearly proves that the much more 
precise corrected quasi-polynomial model is certainly ex- 
tremely accurate, since it is without correction already 
within at most 10~ 2 of the exact result. The next strik- 
ing and important result is the quite remarkable accuracy 
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of the scaling method which is essentially never beyond 
3.10 -3 of the exact result for the MC equation of state 
(the case of the MF-BCS is even better). This gives a 
full validation to the scaling approach. 
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FIG. 3: Reduced square v 2 — ui 2 /ui 2 of the axial mode fre- 
quency in very elongated trap for the Monte-Carlo (MC) 
and the mean-field BCS (MF-BCS) equation of states, as 
a function of k%a, where k% — 1.7N 1 ' 6 /dho is the Fermi 
wavevector of the ideal Fermi gas in a harmonic trap and 
o.ho = (h/mCu) 1 / 2 where u> — (uj x Lj y oj z ) 1 / 3 ' . Experimental re- 
sults are from Ref. [10] . The insert gives the relation between 
k% and fcir(0), used in Fig. 2. 

An essential feature emerging from the results of Fig. 
13 is that the separation between the MC and the MF- 
BCS results is definitely much larger than any theoreti- 
cal uncertainty within the T — hydrodynamic picture. 
It should be consequently possible to distinguish them 
experimentally. In particular the observation of an en- 
hancement of v 2 above 2.5, for kP F a ~ 1, would provide 
evidence for beyond mean-field corrections £4], |25| which 
are absent in MF-BCS. 

Finally we turn to the cigar geometry. We have carried 
out an analysis of the lowest axial compressional mode by 
comparing the results obtained from the scaling method 
with our exact results, in the case studied in Ref. [l7| . 
Again the agreement is perfect. In Fig. [3|we report the 
predictions of the scaling approach applied to the MC 
and the MF-BCS equations of state. In the same figure 
we also show the experimental results from |ldj for the 
axial mode [2(| . It is quite intriguing that these data look 
in better agreement with the mean-field BCS equation of 
state rather than with the Monte-Carlo one. This might 
be due to finite temperature effects. Clearly further ex- 
perimental work would be important to settle this point, 
which is essential for our understanding of the BEC-BCS 
crossover, and in particular to decide which equation of 
state better describes the experimental reality. 
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